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Solar sails employ a unique form of propulsion, gaining momentum from incident and 
reflected photons. However, the momentum transferred by an individual photon is 
extremely small. Consequently, a solar sail must have an extremely large surface area and 
also be extremely light. The flexibility of the sail then must be considered when designing or 
evaluating control laws. In this paper, solar sail flexibility and its influence on control 
effectiveness is considered using idealized two-dimensional models to represent physical 
phenomena rather than a specific design. Differential equations of motion are derived for a 
distributed parameter model of a flexible solar sail idealized as a rotating central hub with 
two opposing flexible booms. This idealization is appropriate for solar sail designs in which 
the vibrational modes of the sail and supporting booms move together allowing the sail mass 
to be distributed along the booms in the idealized model. A reduced analytical model of the 
flexible response is considered. Linear feedback torque control is applied at the central hub. 
Two translational disturbances and a torque disturbance also act at the central hub 
representing the equivalent effect of deflecting sail shape about a reference line. Transient 
simulations explore different control designs and their effectiveness for controlling 
orientation, for reducing flexible motion and for disturbance rejection. A second model also 
is developed as a two-dimensional “pathfinder” model to calculate the effect of solar sail 
shape on the resultant thrust, in-plane force and torque at the hub. The analysis is then 
extended to larger models using the finite element method. The finite element modeling 
approach is verified by comparing results from a two-dimensional finite element model with 
those from the analytical model. The utility of the finite element modeling approach for this 
application is then illustrated through examples based on a full finite element model. 


I. Introduction 

The advent of solar sail spacecraft has brought with it challenges in the areas of materials and manufacturing, in 
the areas of modeling and computation and in the areas of design and control. Much progress has been made during 
the last few years toward completing the vision of solar sail development [1]. Two primary design concepts have 
emerged for a square solar sail supported by four booms radiating from a central hub [2, 3]. These concepts differ in 
the details of sail attachment to the booms, in the design of the booms themselves and in their approach to attitude 
control. One design includes four scalloped-edge sails with multi-point attachments along telescoping booms 
constructed of inflatable isogrid sections. Angular positioning of the sail is accomplished via adjustable vanes at the 
tips of the booms that can be oriented to use the solar pressure to change the solar sail orientation [2]. The other 
design includes sails that are attached to uniform booms at only five locations: the hub and the four boom tips. The 
sails include shear compliant borders to reduce wrinkling [4]. This concept includes a fifth boom and tip mass 
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extending from the central hub out-of-plane with respect to the sail that can be moved to provide attitude control [3]. 
Control simulations with a rigid sail implemented both of these control concepts and established the feasibility of 
these approaches [5, 6, 7]. Recently, the static shape and the dynamic response of a generic sail design similar to 
both concepts have been studied using high-fidelity models and nonlinear finite element (FE) simulations. The effect 
of variation of design parameters on the dynamic response was also considered for both five-point and multi-point 
connection sails [8]. 

The objective of this effort was to develop an initial analytical model to use to consider the impact of flexibility 
on solar sail attitude control and on sail vibration suppression. Ultimately, this question is studied with nonlinear 
transient simulations implementing control strategies on reduced versions of high-fidelity models after the 
simulation methods are verified. The motivational example of a flexible square solar sail design with four booms 
and results from a high-fidelity ABAQUS FE model are given in section II. To enable rapid development of control 
models, a two-dimensional simplified analytical model with only two opposite booms of the four is developed in 
section III. Motion is restricted to one plane. Single-axis torque control is applied to the central hub for vibration 
suppression and attitude control. The distributed parameter equations of motion for this simplified model for first- 
mode motion are derived in section IV. In this paper, several studies are performed with the 2-D simplified model 
whose results are given in section V. First, the effect of disturbance torques with and without linear control feedback 
for harmonic torque distributions that represent sail wrinkling is studied. Second, the controlled response for 
different angular initial conditions is studied. Third, studies are performed to examine the influence of the wrinkling 
disturbance throughout the simulation under no control, delayed control initiation conditions, different feedback 
control schemes, alternate control laws and different angular initial conditions. A solar disturbance model is 
developed and used as input for the 2-D simplified model in section VI, where the solar disturbance model 
determines an equivalent disturbance at the hub attributable to sail shape under solar pressure. Finally, an ANSYS 
FE model is developed in section VII based on the high-fidelity flexible square solar sail FE model of section II. 
Control simulations with the ANSYS FE model are demonstrated in section VIII for unit and random force out-of- 
plane inputs. The simplified analytical model also serves as a verification example for the reduced nonlinear 
transient FE simulations. A summary and directions of continuing work conclude the paper. 


II. Flexible Square Solar Sail 

The analytical model developed for this effort is a simplified representation of a flexible square solar sail. Design 
parameter studies were available in which the variations of the lower modes of a uniform-boom five-point connected 
solar sail were determined for changing design parameters [8]. The effective bending stiffness (El) of each sail 
boom varied from El = 37.2xl0 2 3 to 382xl0 3 Nm as the boom thickness was increased from 4.0 jum to 41 jum. 
Additional information from the parametric study is given within Ref. [8,9]. The first five sail frequencies as a 
function of increasing boom thickness are plotted in Figure la), which is a repeat of Figure 15 of Ref. [8]. The first 
sail frequency as a function of the boom thickness are plotted in Figure lb) for both the five-point and multi-point 
connected sail [9]. For both designs, the shape of the vibration mode associated with the lowest frequency is in the 
shape of a “drum” mode. A drum mode has its maximum deflection at the outer edge of each triangular sail region, 
and the booms move in the same direction as (in phase with) the sail. In essence, the sail and boom motions coincide 
for the first mode. For stiffer booms, the vibration modes associated with the sail membrane become independent of 
the support boom modes. As the development of solar sails continues, it is anticipated that boom stiffness will be 
further reduced, and that the dynamics of these sails will be dominated by combined sail-boom modes and even 
boom-only modes. This observation justifies the study of a reduced flexible solar sail model consisting of uniform 
booms from a central hub with equivalent masses along the boom length to represent the sails. Although not studied 
herein, a non-uniform mass distribution could be considered. 
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Figure 1. a) Sail mode frequency versus boom thickness for a five-point-connected design b) First frequency 
variation for a five- and multi-point connected sail designs versus boom thickness 


III. A Two-Dimensional Analytical Model 

The two-dimensional analytical model is a structure consisting of two opposing cantilevered flexible beams 
connected to a hub. The model includes a central hub mass and two tip masses that are much smaller than the hub 
mass. Figure 2 is a schematic of the model with coordinate definitions. This structure is subject to axial (beam 
longitudinal) and transverse (out-of-plane with respect to the solar sail) accelerations, and v 0 , respectively. A 
control torque, T c , and a disturbance torque, T d , act on the central mass. The mass and rotary inertia of the central 
hub are M 0 and 7 0 , respectively, while M\ and M 2 are the two masses at the tips of two beams. The mass 
distribution, Young’s modulus, and bending stiffness of each beam are denoted £) and J t , respectively, with the 
beam index i = 1,2. The hub inertial rotation is 6 . The axial and transverse displacement functions for the two 
beams are respectively u t (s t , t ) and (s t , t) in which s t is the axial coordinate for beam i and t is the time. 



Figure 2. Schematic of a rigid hub with two flexible beams with end masses including modeling parameters and 

coordinate definitions 


Equations of motion for the flexible 2-D model are developed following examples of a rotating flexible-beam 
satellite in Reference 10 and of multiple-beam flexible solar array support structure in References 11 and 12. 
Standard assumptions are adopted, including the following: 

• the beams are inextensional along the neutral axis, 

• the beams have no shear deformation and warping, 

• the beams are subject to only axial displacement u t and transverse displacement v t , 

• damping is viscous, having the coefficient c t , and 

• gravitational potential energy is negligible. 

From the above assumptions, the axial and transverse displacements are related: 
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where y/f is the bending angle. Third-order Taylor’s expansions for u t and y/i from (1) and (2) are respectively 
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The transverse motion is out-of-plane with respect to the solar sail if we consider this model to represent an 
edge-on view. Following Galerkin’s approach, the motion can be expanded in terms of modes of the underlying 
linear model. Typically, a high-fidelity FE simulation is run for a specific excitation to direct the selection of modes 
for the reduced analytical model. Single mode response can be assumed using the fundamental mode or a higher 
mode, while multiple modes can be assumed to examine nonlinear modal interactions. 

A single-mode model was considered here, based on the first cantilevered bending mode of the boom as shown 
in Figure 3. In addition, the asymmetric, rather than the symmetric, combination of the boom first modes is assumed. 
This is physically consistent with the response expected as a result of an applied torque at the hub. Therefore, the 
transverse (out-of-plane) displacement of the two booms can be expressed in terms of their first bending mode as 

0' = 1 ) and v i( s n t ) = -®\j{si)V l {t) (7 = 2) (5) 



Figure 3. Displacement response actuated by the torque 


IV. Derivation of Equations of Motion via Energy Expressions and Hamilton’s Principle 

The kinetic energy for the model is 

A 

2T = I + Af qI/q + MqVq + J* wi\ (u\ + ) + ^Vq + iq + s^O j ds i + M j [iq ( L j + i/q ] 

0 

A 

+ Ml [vq + iq (a , t)+ J m 2 (zi 2 -Uq ) 2 + (v 0 + v 2 - s 2 of ds 2 

0 

+ M 2 [u 2 (L 2 d)-U()Y + Af 2 [v 0 + v 2 (L 2 ,t)~ L 2 &f 


( 6 ) 


Using Equations 3 and 4 and expanding the kinetic energy in terms of the first mode using Equation 5, the kinetic 
energy is now written as 


T = 


a 0 {u 0 + ci\0 + ci 2 V\ + ^3^i^o cl^Ovq + cl^V^O + d^ViV^UQ + d 2 Vi V\ 


(7) 
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Similarly, the potential energy is 
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According to the extended form of the Hamilton’s principle. 

t sy 2 ^ 2 ^ 


2 Z *2 z 1 r 

SI = S j (T - V )dt + S IJJ yMl - (! 


*' = 1 f, 0 


F 2 *2 

| SW B dt + £ I 


i=l t , 


■ 

J (<2 M , £>’ « I + 6v,.^i )* / 


+ u i ) — v j ^ dt + 


dt = 0 


( 10 ) 


where 8 denotes virtual quantities and are Lagrange multipliers. 5W B is the virtual work associated with motion 
of all the boundaries, and where SW B =(T C + T d )86 in this model. Q u ,Q V are the generalized forces associated 
with the virtual displacements Sui , Svf , where Q u =0, and Q v = -Qv ; - , with Q denoting the damping 
coefficients. 

Inserting Equations 7 and 9 into Equation 10 and considering the virtual work and Lagrange multipliers, two 
nonlinear, coupled equations of motion are obtained: 


1 2 a x 0 + a 5 V x ={T C +T d )-a 4 v 0 
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Equation 11a can be used to solve for the angular acceleration, 

T 
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2 a x 2a x 
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and inserted into Equation 1 lb to yield the nonlinear dynamic equation of motion for the vibration response of the 
idealized sail model 


F, +A,V , + (A 2 - A 3 u 0 )Fj = A 4 V 3 + A 5 V x V 3 + AJT?V X + A 6 v 0 + A 7 (T c + T d ) 


(13) 
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Note that expansion of the motion in terms of the symmetric combination of the first boom modes yields zero 
coefficients for A 3 and A ? , rendering the motion uncontrollable with the resultant control torque T c at the central 
hub. 
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V. Disturbance Control and Angular Position Control Using T c 


Based on the nonlinear dynamic equation of motion (Equation 13), two state variables are 
defined: x x - V\ and x 2 = V\ . First-order state control equations are then written as 


x x = x 2 

< x 2 = y(^3^0 x l +A X 1 3 +^5^l^2 + ^6 X 0 + ^7 _ ^l x 2 - ^2 x l) (14) 

1 — 

Jc = K p x 1 + K d x 2 


The disturbances acting on the hub are T d , Uq , v 0 . These disturbances are the resultants at the hub for excitations 
outside the mission profile. 

A control torque based on the vibrational response states is 

T c =K p x l +K d x 2 t>t x 

in which K p and K d are respectively the control gains for the modal amplitude and modal velocity. From Equation 12 
and with the solution for the vibration rate x 2 , the angular response 6 can be obtained in state form with 

y x = 0 and y 2 = 0 as 


Ti = y 2 

< . T c +T d a 5 . a 4 .. (15) 

y 2 ~ Z — x 2-— x 0 

2 ci x 2,& x 2 & x 

In the following simulations, properties of the two beams are from Reference 12 for a 10-m demonstration solar 
sail or from Reference 8 for a 106-m boom solar sail. Tables 1 and 2 summarize the properties and resulting natural 
frequencies for the first three modes of the beams. 

Table 1. Single-Boom Properties for the 10-Meter Analytical Model 


Properties 

Value 

Bending stiffness (El) 

2660 Nm 2 

Mass per meter 

70 g/m 

Length 

1 m 

Hub mass 

2.5 kg 

Tip mass 

0.1 kg 

Damping coefficient 

0.015 

First three natural frequencies 

1.639, 28.448, 51.389 Hz 


Current interest in the effects of wrinkling [14] lead to questions about effects of responding sail shape, 
prompting modeling of the sail disturbance as a harmonic function resulting from the shape of the sail as it responds 
in the first mode as 
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T d =0.1cos(10.30 t<t e and w 0 = 0 v 0 = 0 


For now, axial and thrust components of such a disturbance are assumed to be small in comparison to the torque 
disturbance, and are neglected. Future efforts will combine time-varying disturbances resulting from sail shape (see 
the next section for preliminary study) with this simulation. 

Table 2. Single-Boom Properties for the Integrated Sail Model 


Properties 

Value 

Modulus 

1.24x10" N/m 2 

Density 

1908 kg/m 3 

Poisson’s ratio 

0.3 

Outer radius 

0.22862 m 

Wall thickness 

7.5xl0~ 6 m 

Length 

106 m 

Hub mass 

2.5 kg 

Tip mass 

0.1kg 

Damping coefficient 

0.005 

First three natural frequencies 

0.060, 0.378, 1.066 Hz 


The first MATLAB [15] control simulation focused on the modal amplitude response induced by a disturbance 
torque without and with a feedback control torque. The simulation uses the properties from a 10-m demonstration 

solar sail. For this simulation, the initial conditions are Vi(0) = 0m, f\(0) = 0. m/s, ^(0) = 0 rad, and 
6^(0) = Orad / s , the disturbance torque is T d = 0.1cos(10.3/) where t<l.5s and the control torque is 
T c = K d V ( t ) , where t>0s and K d = 1 N -s . 

Figures 4 and 5 present the results for two simulations. The first is uncontrolled, with the disturbance ending 
after 1.5 seconds. The second simulation has the same disturbance, but also includes a modal velocity feedback 
control applied for the entire simulation. The control is seen to damp the transient vibrational (modal amplitude) 
response at the beginning of the simulation and again when the harmonic disturbance discontinues. With larger 
gains, the vibration suppression is more immediate, as expected. Figure 6 is the same simulation, but with the gain 
eight times larger. 
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Figure 4. Modal amplitude response with disturbance T d = 0.1cos(10.3f) where t <1. 5s and no feedback control, 
x 10' 4 



Figure 5. Modal amplitude response with disturbance T d - 0.1cos(10.3£) where t<l.5s and 
feedback control T c = K d V(t ) , where t > Os and K d =IN -s . 
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x 10' 4 



Figure 6. Modal amplitude response with disturbance T d - 0.1cos(10.3£) where t < 1.5s and feedback control 

T c = K d V(t ) , where t > Os and K d - 8N • s . 

The second simulation focuses on the coupled angular and modal responses resulting from an angular rate initial 
condition. This simulation is without solar radiation or other disturbance torques. The example uses properties from 

a 106-m boom integrated solar sail. For this simulation, the initial conditions are Fj (0) = 0m , Fj(0) = 0m/ s , 

Q\ (0) = Orad , and ^(0) = Ora d/s , and the control torque is T c - K a (Q(t) - 0$) + K b @(t), with K a =-1000, 

K b = -100 . The control law is a linear feedback control to move the sail to the reference position (0 = 0) and 
reduce the angular rate to zero. Results are presented in Figures 7 through 9. 
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Figure 7. Controlled response of a 106-m boom model from angular rate initial conditions: angular response. 
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Figure 8. Controlled response of a 106-m boom model from angular rate initial conditions: modal response. 



Figure 9. Controlled response of a 106-m boom model from angular rate initial conditions: control torque. 

The next simulations change the scenario: the disturbance now persists throughout the 10-second simulation and 
the control is initiated at time t = 1.5 seconds. Furthermore, the model damping is decreased to 
c \ = c 2 = 0.0057V -s /m . Figures 10 and 1 1 present the modal amplitude response and the angular position response, 
respectively, for no control. 

Figures 12 and 13 present the modal amplitude and the angular position with a feedback control with only modal 
amplitude (no velocity) feedback. Note that to have any influence, the gain for the position feedback (modal 
amplitude) must be large since the modal amplitudes are very small. To further illustrate this point, the same 
simulation is run, but with the addition of a small velocity feedback term. Figures 14 and 15 are plots of the modal 
amplitude response and the angular position response, respectively. A significant reduction of the modal amplitude 
is seen in Figure 14 compared to Figure 12. The results of both controlled simulations show a linearly increasing 
angular position after the control is initiated. 
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To control the angular position of the flexible solar sail, a different control law is implemented. For these 
simulations, the harmonic disturbance is not present. First, the orientation angle is controlled to move the solar sail 
from an angular position of 0° to an angular position of 5°. The control torque is a combination of a constant torque 
to initiate the motion and an angular rate feedback to stop the maneuver. Figures 16 and 17 present the modal 
angular position response and modal amplitude time histories for this maneuver, respectively. Their control torques 
are represented as below 

0 t < 2s, 


< 0.2 Nm t > 2 s and 6 < 4.03° 

K a 0(t) 0 > 4.03° ,K a = -50 Nms/rad 


In the final simulation of this section, the angle 6$ is controlled in a maneuver from 8° to 0° to -8° with control 
torque designed as a proportional feedback of the angular position error and the angular rate, 
T c - K a ( 6(t ) - #o) + K b 0{t) with K a = -20Nm / rad , K b - -50Nms / rad . Figures 1 8 through 20 present 
respectively the angular position, the boom tip mass response and the torque time histories for this maneuver. Note 
that if the gain K a is less than zero and greater than -20 Nm / rad , the control torque would be smaller, but the 
control time would be longer. 


x 10' 4 



Figure 10. Modal amplitude response with disturbance T d - 0.1cos(10.3£) where t < 10s and no control. 
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x 10' 3 



Figure 1 1 . Angular position response with disturbance T d =0.1cos(10.3 1) where ^<10^ and no control. 


x 10' 4 



Figure 12. Modal amplitude response with disturbance T d = 0.1cos(10.3£) , where t < lOs and feedback control 
T c = K p V(t) + K d V(t ) , where t > 1.5s , K p = 60000A , K d = ON -s . 
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Figure 13. Angular position response with disturbance T d - 0.1cos(10.3/) , where t < 10s and feedback control 
T c = K p V(t ) + K d V(t ) , where t>l.5s, K p = 60000 A , K d = ON -s. 


x 10' 4 



Figure 14. Modal amplitude response with disturbance T d - 0.1cos(10.3£) , where t < 10s and feedback control 
T c = K p V(t) + K d V(t ) , where t > 1.5s , K p = 60000 A , K d =2N-s. 
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Figure 15. Angular position response with disturbance T d - 0.1cos(10.3/) , where t < lO.s and feedback control 
T c = K p V(t ) + K d V(t ) , where t>l.5s, K p = 60000TV , K d = 2N • j . 



Figure 16. Angular position response for angular positioning maneuver from 0° to 5° with control. 
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Figure 17. Modal amplitude response for angular positioning maneuver 6 0 from 0° to 5° with control. 



Figure 18. Angular position response for angular positioning maneuver 9 {) from 8° to 0° to -8° with control 
T c - K a ( 6{t ) -0 O ) + K b 6(t) where K a - -20 Nm / rad , K b - -50 Nms / rad . 
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x 10' 3 



Figure 19. Boom tip displacement response for angular positioning maneuver # 0 from 8° to 0° to -8° with control 
T c = K a {0(f) - <9 0 ) + K b 6(t) where K a = -20 Am / rad, ^ = -50Ams / rad . 



Figure 20. Control torque for angular positioning maneuver from 8° to 0° to -8° with control 
T c = K a (0(t) - 0 O ) + where K a = -20 Am / rad, K b = -50 Nms / rad . 
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VI. Solar Radiation Disturbance 

Disturbances acting on the hub can be specified independent of the sail response or can be computed based on 
the response shape of the solar sail. An algorithm was developed and implemented in MATLAB to determine the 
solar radiation pressure disturbance effects. The algorithm first divides the phenomena model sail into n elements 
each of equal length M, and then calculates the displacement of each element from the initial flat sail plane ( D ) 
according to the mode shape function corresponding to the current time step. The distance from the center of each 
element a to the center of the sun (r ) is calculated using the following formula: 

r(d) = R- [(Af ((0.5» + 1) - a) - 0.5M]cos(S) + 5(a) (16) 

where the constants are as follows: R = 150xl0 9 meters is the distance from the Sun at 1 AU, where the solar 
intensity is /= 4x1 0 26 watts, and the speed of light is C = 3x1 0 8 meters per second, S is the counter-clockwise angle 
from sail plane to line connecting the sail hub to the center of the Sun, and W is the width of each element in meters. 
See Figure 21 for geometric definitions. 

The MATLAB script calculates the solar intensity (/ in N/(m-s)) and force of the Sun (f in N/m 2 ) due to that 
intensity at each element a using the following equations 


«■(«) = — m= 7 ^ (17 a,b) 

4 nr{a) c 



As shown in Figure 22, the Sun force at each element is calculated by the projection of each element onto the 
sail plane normal to the Sun. The analytical model shape is approximated by a discrete number of straight line 
members. The slope of each member is used in the torque calculation, but not in the force calculation below to 
ensure that any disturbance torque results from the differential in the distance of opposite parts of the sail to the Sun, 
and not the approximation. The Sun force (F) is therefore calculated in Newtons from Newtons per square meter in / 
above. The slope of each element is used to calculate the angle from the element normal to the sail plane ( eleangle ), 
where 
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(18 a, b) 


F(a) = f{a)m sin(S)fF 

eleangleia ) -n- ^ - a.YCta,n(eleslope(a)) 


The perpendicular moment arm of the force vector acting along the element normal is calculated by first finding 
6 , the angle from the sail plane normal to the line from the center of the element to the center of the sail, £ , where X 
is the x-distance from the center of the sail to the center of the element. The angle P is the angle from £ to the 
“bottom” of the right triangle, or the line connecting the perpendicular moment arm to the center of the element, 
where 


0 - arctan 



if eleangle(a ) > k/2 

P (a) = n - atan (D(a)/X(a)) - eleangle(a ) 
elseif eleangle(a ) < k/2 

P (a) = atan (D(a)/X(a)) + eleangle(a ) 
else 

P (a) = tt/ 2; 
end 


The length of the perpendicular moment arm is calculated as sin(/?) multiplied by the length of £ . To enforce 
the correct sign on the moment calculation, the length of the moment arm is multiplied by the sign of X in Equation 
19. This ensures that forces on opposite sides of the sail hub produce opposite signed torques. 

elemomentcrm{a) = ^X(a) 2 +D(af sin {J3{a)) sign(X(a)) (19) 

The total disturbance torque T d is then the sum of the force on each element F is multiplied by the perpendicular 
moment arm summed over the entire sail. 

n 

T d =y] F(a) elemomentarm(a ) (20) 

a = 1 



Figure 22. Solar radiation disturbance computation parameters 
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A study of solar radiation effects was conducted first to examine the orientation of the sail with respect to the 
sun, the amplitude of the modal response and trends in the resulting disturbance torque. The geometric and material 
properties used for this study were for a 10-m test article sail, allowing efficient computation of the trends in 
contrast to the time required to compute trends with a larger sail model. The trends seen here are valid for larger 
models. Figure 23 presents the results of this study, and we see that the largest disturbance torques result when the 
sail plane is oriented 30 degrees from perpendicular to the sun line. When the sail plane is perpendicular to the sun 
line (at 0 degrees), the disturbance increases with the amplitude of the modal response. When the sail is oriented 
parallel to the sun line (at 90 degrees), the disturbance effects are independent of the shape of the sail response. 



Figure 23. Disturbance torque at different angles between solar sail and sun vector 

(7-m boom) 


The next simulation focused on the angular response induced by solar radiation disturbance torque resulting from an 
initial modal displacement condition. The example uses properties from a 10-m demonstration solar sail with a 
boom length oil m. Due to the order of magnitude smaller size of the sail in this model, the radiation intensity is 
increased accordingly to 4xl0 36 w. The distance from the Sun is 1.5xl0 n m; the speed of light is 3xl0 8 m/s. For this 

simulation, the initial conditions are 1^(0) = 0.1m, 1^(0) = 0m/s, 6i(0) = 0 rad, and 6^(0) = Orad / s , the 
disturbance torque is computed using the solar radiation disturbance model and the control torque is 
T c - K a (0(t) - Oq) + K b 6(t) with K a =-100 ,K b =-10. The angular response without control is seen in Figure 
24. After a long time, a change of angular orientation is seen. The angular response with control is seen to diminish 
in Figure 25a), along with the control torque in Figure 25b). 
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0.025 


0.02 



Figure 24. Uncontrolled response of a 10-m sail to solar radiation disturbance. 




b) Control torque 

Figure 25. Controlled response of a 10-m sail to solar radiation disturbance. 
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The next simulation presented focused on the angular response induced by solar radiation disturbance torque 
resulting from an angular rate initial condition. The example uses properties from a 106-m boom solar sail. The 
initial angle between the sun and solar sail surface is 30 degrees. For this simulation, the initial conditions are 

^(0) = 0 m , V\(0) = 0ml s , #i(0) = 0 rad , and ^(0) = 0.1 rad Is , the disturbance torque is computed using the 
solar radiation model and there is no control torque applied. Figure 26a) presents the boom tip displacement of the 
uncontrolled simulation, while Figure 26b) presents the angular displacement. Note that the steady state angular 
orientation is not zero. 



a) Boom tip displacement 



b) Angular position 

Figure 26. Uncontrolled response of a 106-m boom sail to solar radiation disturbance. 


The next simulation is the same as the previous one, but with control. It focused on angular response induced by 
solar radiation disturbance torque resulting from an angular rate initial condition. The example uses properties from 
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a 106-m boom solar sail. The distance from the Sun is 1.5x10 n m; and the speed of light is 3xl0 8 m/s. For this 
simulation, the initial conditions are ^(0)^0 m, F 1 (0) = 0 m/s, &i(0) = 0 rad, and ^(0) = O.lrad / s , the 
disturbance torque is computed using the solar radiation model and the control torque is 
T c - K a (0(t) - 0) + K b 6{t ) , where K a = -100 ,K b = -100 . Figure 27 presents the results of this simulation. 


0.1 . T . T . T . T . 
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0.06 ■ 
e; 
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□ 
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0 

. 0.02 ' ' * ‘ * • * * * 
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a) Boom tip displacement 



b) Angular position 



c) Control torque 

Figure 27. Controlled response of a 106-m boom sail with solar radiation disturbance: modal, angular, and control 

torque. 


VII. Finite Element Modeling Approach 

The complicated dynamics of large solar sails may best be modeled using the finite element method. However, 
use of high-fidelity FE models in transient controls simulations can be computationally intensive. A study involving 
use of reduced FE models in attitude control simulations for the solar sail application has been undertaken. The 
method was implemented using the ANSYS [16] FE software and the MATLAB numerical analysis software 
package. The FE software is used to generate a structural system model in standard second-order matrix form. The 
system model is then reduced via Guy an reduction [17]. The reduced model is then imported to MATLAB. A 
representation of a control system is incorporated with the system model, and the model is converted to first order 
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form, and simulated using built-in MATLAB numerical integration functions. The overall approach is similar to the 
approach outlined in Reference 18 for studying vibration control of turbomachinery blading. 

The second order system FE model, including the control system and possible external disturbance loads, can be 
written as 


[M] {d} + [C] {d} +[K] {d} = {F e } + {F c } (21) 

where {d} is a vector of nodal displacements, {j}is a vector of nodal velocities, {j}is a vector of nodal 
accelerations, {F e } is a vector of external disturbance loads, and {F c } is a vector of control forces and moments. 
\M\, [C], and [K] are the structural mass, damping, and stiffness matrices, respectively. 

The model can be reduced via Guyan reduction to a set of m master degrees of freedom (DOF), and written, 

[M r ]{d m } + [C r ]{d m } + [K r ]{d m } = {F e _ r } + {F c _r } (22) 


where {d m } is a m x 1 vector of nodal displacements of the master DOF, { d m } is a m x 1 vector of nodal velocities 

of the master DOF, {d m } is m x 1 vector of nodal accelerations of the master DOF, {F e . r } is a m x 1 reduced vector 
of external disturbance loads, and {F c . r } is a m x 1 vector of control forces and moments. Matrices [M r ], [CJ, and 
[K r \ are the reduced structural mass, damping, and stiffness matrices, respectively, each of size mxm. 


To convert the reduced model to state-space form, a vector of state variables, {q}, is defined, where 


{<?}= 



The following matrices are also defined: 

[Z] [/] 

[K r ] [C r ] 


[*] = 


-[/] [ Z ] 

[Z] [M r \ 


[H]= 




(23) 


(24) 


where, for a model with m retained master DOF, [Z] is an m x m matrix of zeros, [7] is an m x m identity matrix, and 
{Zj is an m x 1 column vector of zeroes. Using Equations 22, 23 and 24, the model can be written 

[R]{q}+[H]{q} = {F e _ r }+ {F c _ r } (25) 

and then rearranged into standard state-space form as 

{q} = [A]{q} + [B]{u} + {P e } = -[R]- l [H]{q} + [R]- l {F c _ r } + [R]- l {F e _ r } (26) 

The vector, {. P e }, relates to external disturbance loads. In Equation 26, for the systems studied here, the 
coefficient matrix, [A], is of size 2m x 2m. In general, for a multi-input system, {u} is a vector of control inputs. In 
typical feedback control applications, the control input is proportional to measured system responses. It is assumed 
here that the control input, {«}, can be written as 

{u}=-[G]{q} (27) 

where [G] is a matrix of control gains. Substituting Equation 27 into Equation 26 results in the closed-loop system 
equations, 

fe}=[[^]-[5][G]]fe}+{F e } (28a) 


or, 

\Acl ]M + { p e } (29b) 

The system response, { q(t )}, due to loads in the vector, { P e }, and/or some initial conditions, {g(0)}, can then be 
calculated by standard numerical methods using this set of system equations. In this work, the built-in MATLAB 
function for performing Runge-Kutta integration is used. 
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VIII. Finite Element Simulations 


In an initial study, to gain preliminary verification of these FE modeling and simulation procedures for the solar 
sail application, a FE model was created corresponding as closely as possible to the two beam distributed parameter 
model described above. The model consists of twenty ANSYS 2D beam elements (ANSYS Beam3 elements), each 
with two end nodes and three DOF per node. The DOF at each node are translations, U x and U y and rotation, 6 Z . 

Three lumped-mass/inertia elements (ANSYS Mass21 elements) are included at the hub and boom ends. The hub is 
constrained to zero translation. The reduced model retains as master DOF the rotation at the hub, 6 Z , and the 
translation, U y , and rotation, 0 Z , at the boom ends. This example is a 10-m demonstration solar sail. A schematic 

of the model and required modeling parameters are shown in Figure 28. For this simulation, an initial rigid body 
angular velocity is applied to the structure of 1.0 rad/s (counter-clockwise, according to Figure 28). In order to 
implement an initial rigid body rotation for the structure, initial angular velocities are imposed at the retained master 
6 z DOF at the left end, right end, and hub. Also, corresponding initial translational velocities are imposed at the 
retained master U y DOF at the left and right ends. So, for an initial rigid body angular velocity of 1 rad/s (counter- 
clockwise), there is an initial translational y-direction velocity of -7 m/s at the left end, and +7 m/s at the right end. 
There is no disturbance torque assumed, and the control torque is defined as 

T c =K a {0{t)-e o ) + K b e{t) where K a =-1000 ,K b =-100. 

Figures 29 and 30 present the results of these simulations. When the uncontrolled and controlled responses for 
the 2-D analytical model in Figure 29b) are compared to those of the reduced FE model in Figure 30b), a slightly 
different response is observed. The reduced model response is more lightly damped, and in particular, the trough 
following the initial peak is deeper and delayed in Figure 30b). This is not surprising as more modes participate in 
the reduced finite element model simulation than in the 2-D analytical model. However, the results from both 
models are in good agreement, which provides some verification for the finite element modeling procedures. 

Additional simulations for the reduced FE modeling approach were conducted based on a model constructed to 
be of a size more typical of that expected for the integrated solar sail models. This model was derived from the high- 
fidelity model used to generate Figure 1 and is described in References 8 and 9. The FE mesh is seen in Figure 31, 
where node 1 is at the hub, and nodes 2, 52, 102, and 152 are at the tips. Details of the model are as follows: 

• Four Booms : properties given in Table 2; ANSYS Pipe 16 elements - 50 elements per boom. 

• Control Mast : Length = 2 m (directed out of plane of sail from hub); 10 ANSYS Pipe 16 elements; radius = 

0. 01705 m ; wall thickness = 0.005 m ; modulus of elasticity = 1.24xlO n TV/m 2 ; density = 7660 kg/m 3 . 

• Membrane : ANSYS Shell41 elements (connected along booms); shell thickness = 2.5 xlO" 6 m; modulus of 
elasticity = 2.48x1 0 9 N/m 2 ; density = 1572.45 kg/m 3 . 

• Lumped Masses : Spacecraft - 228 kg at end of control mast; stowage - 54.84 kg at hub; tip masses at ends 
of booms - 0.58 kg. 

• 18 Retained Master DOF : At hub: U x ,U y ,U Z ,0 x ,6 y ,0 Z ; at each of the four boom tips: U x ,U y ,U Z . 

• Structural Damping : Proportional to the stiffness matrix, [ K \ , so that the damping matrix, [C], is defined as 
[C\=0.001*[K\. 

• Control Moments : Applied at the hub, and based on feedback of hub angular displacements and angular 
velocities, according to: M x =-50,000 0 X - 5,0000^,717^ = - 50,000 0^ - 5,0000^ . 

M x is a control moment at node 1 about the x-direction, M y is a control moment at node 1 about the y- 
direction, 6 X and 6 X are the angular displacement and velocity, respectively, about the x-direction at node 

1, and 0 V and 6 y are the angular displacement and velocity, respectively, about the y-direction at node 1 . 

In an initial simulation using this model, the excitation was a unit force applied for 0 < t < 1 s in the z-direction 
(out-of-plane) at two adjacent boom tips (nodes 2 and 52). The system response was simulated for the controlled and 
uncontrolled case. The results for rotation about the x-direction at the hub (node 1) are shown in Figure 32. In these 
results, linearly increasing rotation is seen in the uncontrolled response, with flexible solar sail motion seemingly 
superimposed on this underlying response. Hence, control designs should consider both rigid and flexible motion of 
the solar sails. 
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A second simulation was run with the same model, but with a disturbance due to random out-of-plane forces. 
Disturbance force magnitudes ranged from 0 to IN. Different forces were applied at each boom tip (nodes 2, 52, 
102, and 152) and the forces were updated every 0.001 seconds. The random input was implemented using the 
MATLAB “rand” command. Figure 33a) presents the uncontrolled and controlled results for angular displacement 
about the y-direction at the hub, and Figure 33b) presents the angular displacement about the x-direction at the hub. 
Simulations with three-dimensional solar radiation disturbance models integrated with the reduced finite element 
model are possible using this approach, but were considered beyond the scope of this effort. 


Modulus, E=3,192,000 N/m 2 
Density, p = 0.7 kg/m 3 



Figure 28. Finite-element model for verification of the reduced model approach. 




a) Uncontrolled angular response b) Controlled angular response 

Figure 29. Analytical model results. 
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Angular Displacement (rad) 




a) Uncontrolled angular response 


b) Controlled angular response 


Figure 30. Reduced FE model results. 



Figure 31. Representative finite-element model of a solar sail. 
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Figure 32. Simulation results with representative finite element solar sail model. 


9 y at Hub with Random OP Force Excitation at Boom Tips 



^ -3 6„ at Hub with Random OP Force Excitation at Boom Tips 
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a) Angular response 6 y b) Angular response O x 

Figure33. Simulation results for a random force excitation at the boom tips with representative solar sail finite 

element model. 


IX. Summary and Future Work 

In this paper, three models were developed for the response of flexible solar sails. The first model is a two- 
dimensional idealization of a flexible solar sail representative of the class of solar sail systems in which the sail and 
booms move synchronously for the lower modes. The model consists of two flexible booms with a central hub and 
tip masses. Equations of motion were developed for the coupled vibrational response (reduced to the response of the 
first asymmetric mode only) and angular motion. Control simulations were run with linear feedback control of the 
vibrational response states and of the angular response states to reduce vibrational motion and to execute 
repositioning maneuvers. This model serves to develop physical understanding of the influence of flexibility. In 
addition, this model is used to validate planned closed-loop control simulation implementations with commercial 
finite element software. The second model is a two-dimensional solar radiation disturbance model to calculate the 
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effect of solar sail shape on the resultant thrust, in-plane force and torque at the hub. The torque from this model was 
then input to the idealized model to demonstrate vibration and attitude control. A reduced finite element model also 
was developed and used to demonstrate vibration and attitude control of the solar sail. Limited verification of the 
response of the two-dimensional finite element model was provided by comparison to the response of the first 
analytical model developed herein. The finite element modeling approach was then extended to simulate both 
controlled and uncontrolled response of a full finite element model of a solar sail. 
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